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Relations between connected and self-avoiding walks in a 

digraph 

Thibault Espinasse* and Paul Rochet^ 


Abstract 

Walks in a directed graph can be given a partially ordered structure that extends to 
possibly unconnected objects, called hikes. Studying the incidence algebra on this poset 
reveals unsuspected relations between walks and self-avoiding hikes. These relations are 
derived by considering truncated versions of the characteristic polynomial of the weighted 
adjacency matrix, resulting in a collection of matrices whose entries enumerate the self¬ 
avoiding hikes of length £ from one vertex to another. 
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1 Introduction 

A directed graph G = (V, E) is defined by a vertex set V = {ui,..., uat} and set E of ordered 
pairs of vertices representing the directed edges of G. Directed graphs, or digraphs, have been 
extensively used in the literature as mathematical models to describe actual phenomena such as 
social interactions [5], road traffic |9], physical processes jl] or random walks |3] among many 
others. In most models, it is convenient to allocate weights to each edge of the graph in order 
to incorporate some additional information. For example, a weight ujij between two nodes Vi,Vj 
may serve defining the transition probability of a random walk, the speed limit or the type of 
road in the traffic network and so on. A finite graph on N vertices is then characterized by 
its weighted adjacency matrix W = which accounts for the level of interaction 

between vertices. In this way, every edge of G is identified with a weight uJij. 


An oriented walk on the digraph G is defined as a succession of contiguous edges uiij. Seeing 
the weights uiij as formal variables, a walk w of length i from vi to Vj can be viewed as a degree 
i monomial w = The weighted adjacency matrix then provides a practical tool 

to handle walks on the graph, as they can be derived from analytical transformations of W. For 
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instance, the (i,j) entry of W^, given by enumerates all walks of length 2 

from Vi to Vj. The introduction of the weighted adjacency matrix W to describe the walks on a 
graph goes back to the 60’s. In |10) and the spectral properties of a graph are investigated 
via the determinant and characteristic polynomial of W. Digraphs also provide a useful tool to 
compute the determinant and minors of sparse matrices, as discussed in |15| . For general results 
on spectral graph theory, we refer to laiH]. 

In |17) . the author shows that the coefficients of the characteristic polynomial of W can be 
interpreted in term of the self-avoiding cycles in G. In this paper, we derive a similar result con¬ 
cerning self-avoiding hikes, defined as a generalization of walks to possibly unconnected sequences 
of edges. We construct a collection of polynomials of W whose entries enumerate the self-avoiding 
hikes of a given length from one vertex to another. The polynomials are obtained as Cauchy 
products of the characteristic polynomial coefficients with the sequence of successive powers of W. 

The analytical expression linking the self-avoiding hikes and the walks on the graph hides a 
deeper connection when considering each hike individually. Precisely, the relation can be investi¬ 
gated in the partially ordered set formed by the hikes. In this context, combinatorial properties 
arise when studying functions of the hikes in the reduced incidence algebra of this poset. In 
particular, we show that the number of different ways to travel a closed hike can be expressed 
in term of its self-avoiding divisors via a Mobius-like inversion on this poset. Another result on 
the decomposition of a walk into self-avoiding components is then derived. 

The paper is organized as follows. Definitions are introduced in Section as well as the 
preliminary result. Section is devoted to the study of the different relations between self¬ 
avoiding hikes and walks, where many combinatorial properties are investigated. The results are 
verified on specific examples in Section 

2 Notations and preliminary results 

Let G = {V, E) be a labeled directed graph, or digraph, with finite vertex set V = {ui,..., vn} 
and edge set E which may contain loops. The adjacency matrix of G is defined as the N x N 
matrix A with entries aij equal to one if Vi is connected to Vj and zero otherwise. Because G is 
directed, A may not be symmetric {vi can be connected to Vj without Vj being connected to Vi). 
In this paper, an edge always refers to a directed edge. 

The adjacency matrix can be used to derive numerous properties of a graph. For instance, 
the (i, j) entry of A^ gives the number of walks of length £ from Vi to Vj. When one is interested 
in each walk specifically, a useful tool is to allocate a weight, or variable, to each non-zero entry 
of A. In this way, the digraph G is characterized by its weighted adjacency matrix 

where the Wjj’s are real variables, setting uiij = 0 whenever there is no edge from Vi to vj. An 
edge of G can then be identified with a non-zero variable ujij and a walk w from Vi to Vj with the 
product w = of the edges composing it (two walks are thus considered equal if 
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they are composed of the same edges, counted with multiplicity, regardless of their order). The 
walk w is closed (a cycle) i = j and open (a path) otherwise. Moreover, w is simple if it does 
not cross the same vertex twice, that is, if the indices i,ii, j are mutually different, with 

the possible exception i = j if w is closed. Loops uiu and backtracks UijOOji are considered cycles 

of length 1 and 2 respectively. 

The representation as monomials in the formal variables ujij provides a simple multiplicative 
structure on walks. The cycle-erasing procedure of Lawler |13) shows that a walk from Vi to Vj 
can always be decomposed as the product of a simple walk from Vi to Vj and cycles, as illustrated 
in Figure However, the reverse is not true in general as the product of a simple walk and 
simple cycles might not be connected. In this paper, we define a new object, called hike, which 
extends the definition of a walk by relaxing the connectedness condition. 

Definition 2.1 A hike from Vi to vj is a monomial h = that can be decomposed into 

the product of a simple walk from Vi to vj and simple cycles. 

Properties of walks naturally extend to hikes. A hike is closed if it is a product of simple 
cycles and open otherwise. By convention, the trivial cycle 1 is considered a closed hike of length 
0. Similarly as for walks, a hike is self-avoiding if it does not cross the same vertex twice. The 
connected components of a self-avoiding hike are simple and vertex-disjoint. In this setting, a 
walk can be viewed as a connected hike. 

In the sequel, the length of a hike h (its degree) will be denoted by i{h) and its number of 
connected components by n{h). A walk is a connected hike, that is, a hike h such that n{h) = 1. 
Moreover, we denote by V{h) the set of vertices crossed by h and |H(/i)| its cardinal. Remark 
that an open hike h is self-avoiding if, and only if, |R(/i)| = £{h) + 1 while for a closed hike, 
self-avoiding is equivalent to |R(/i)| = £{h). 



Figure 1: The walk w = uji 2 U! 23 UJ 34 U! 42 W 25 from vi Figure 2: The self-avoiding hike h = UJ 12 OJ 34 UJ 45 UJ 33 
to V 5 is the product of the simple walk ijJi 2 UJ 25 and from vi to V 2 contains n{h) = 2 connected compo- 
simple cycle uj 23 W 34 UJ 42 . nents: the path uji 2 and the cycle UJ34UJ45UJ53. 

The following notations are used throughout the paper. 

- W denotes the set of walks. 

- H denotes the set of hikes (open and closed). 
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- C denotes the set of closed hikes. 


- S denotes the set of self-avoiding hikes. 

For these sets, we may specify the end-vertices in index and/or their length in exponent, e.g. '}i\^ 
is the set of hikes of length i from Vi to vj. Remark that W, C and S are snbsets of Ti, while Su 
and Wii are snbsets of C for all i = 1, ..., N. The set of self-avoiding closed hikes is written as Cn5. 

It is known that the t'-th power of W ennmerates with mnltiplicity the walks of length i 
on the graph. The (i, j) entry of is a homogenons degree £ polynomial attribnting to each 
walk w of length i, a coefficient fij{w) connting the nnmber of ways to write w as a snccession 
of contignons edges starting from Vi and ending at Vj (the fnnction fij is compnted on some 
examples in Fignrej^. The formal series associated to the fnnctions fij follows by the identity 

(I-W)-.i = (l + W + W2 + ....)^^.= Y.Mw)w, (1) 

loSW 

which holds whenever |||W||| := sup^gj^jv ||^||=i ||Wn|| < 1, with ||.|| the Enclidean norm in 
Remark that if w is an open walk, there is at most one conple (i, j) for which fij{w) is non-zero. 
This property is no longer verihed for a closed walk c in which case fn{c) may take different 
positive valnes for different nodes Vi G V{c), as illnstrated in Fignrej^ 


@ 


© 
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Fignre 3: The walk w = ujiqU!q3U!32(^29<-^95(^54(^49(^97<-^76^69^98 (gray and black edges) has fis{w) = 6 
contiguous representations starting from vi and ending in ug. Each representation relies on an ordering 
of the simple cycles a;g3W32a;29, Wg5a;54u;49 and ujgTUtTQUjQg to travel w. The cycle c = tLig 5 a; 54 W 4 ga;g 7 a; 76 W 69 
(black edges) has f 99(c) = 2 contiguous representations starting from vg and fii{c) = 1 contiguous 
representation for any other vertex vt G T(c) \ {ug}. 

Closed hikes also play an essential part in analytical graph theory. A simple reason is that a 
self-avoiding closed hike can be associated to a permntation cr on a snbset of {1,..., Af} by the 
relation Ca = Using that the nnmber of connected components n(co-) is linked 

to the signatnre of a throngh the identity sgn(cr) = (—1)'^“"'^'^), we obtain an expression of the 
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determinant of W by 

det(W) = ^sgn(cj) = (-1)^ Y1 

ceCricS^ 

where the first sum is taken over all permutations over {1,A^} and the second sum over all 
self-avoiding closed hikes of length N. A more general formula, given in Theorem 1 in links 
the coefficients of the characteristic polynomial 


N 

x(A) = det(AI-W) = 

k=0 

with the self-avoiding closed hikes of length k by 

^^0 = 1, i’k= Yl k = l,2,... (2) 

ceCn^'" 

the coefficient V’fc being trivially zero for k > N. 


The coefficient ( —is reminiscent of a Mobius function. In fact, the function ^ defined 
over C by //(!) = 1 and 

if g ig self-avoding 
0 otherwise. 


fi{c) := 


( 3 ) 


is the Mobius function of the trace monoid of a partially commutative version of closed hikes, 
studied in |6] under the name circuit. The particular value of the characteristic polynomial 


N 

x(l) = det(I -W) = Y'^k = Y 


k=o cec 


( 4 ) 


will turn out to be particularly important for our purposes, as it gives the formal series associated 
to /i. Although the definition of [i is restricted to closed hikes, some of its properties have direct 
repercussions on open hikes. This is due to the fact that an open hike h between two different 
vertices Vi, Vj can be expressed as a closed hike to which the edge ujji has been removed. Actually, 
one verifies easily the following equivalence for i ^ j 


h G kiij \ y hujji G C (5) 

where we recall that %ij is the set of hikes from Vi to vj. This property provides a natural 
extension of /r to open hikes. For all directed pairs of vertices (i,j), define 

( (—if i ^ j and h G Sij 

IJ.ij{h) := —nihojji) = < (—if i = j, h G C n 5 and Vi ^ I^(h), (6) 

I 0 otherwise. 


If i ^ j, the function mj only takes non-zero values for self-avoiding hikes from Vi to Vj. 
In particular, = 0 and ^ij{h) = 1 if /i is a simple path. On the other hand, is 
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non-zero only if h is closed and does not cross Vi. We have for instance = ^i{uii) = 1 (the 

convention = 1 instead of —1 jnstihes the minns sign in the dehnition of ^ij). Remark 

moreover that fJ,ij{h) is nnll for all hikes h of length i{h) > N since a self-avoiding hike on the 
graph has maximal length N. 


The main resnlt of this paper (Theorem 3.2) exhibits a dnality between the fnnctions /, 


iij 


connting the nnmber of connected representations of a hike and the fnnctions Hij whose snpports 


are restricted to self-avoiding hikes. The dnality is actnally a conseqnence of Lemma 2.2 which 
expresses the matrix formal series M = of fiij, dehned by 


rriin : = 


i,j = 1, ...,iV, 


(7) 




in fnnction of the powers of W. 


Lemma 2.2 //|||W||| < 1, 

l 

( 8 ) 

£>0 k =0 

where the V'fc, k = 0,1, ... are the coefficients of the characteristic polynomial defined in ([^. 

Proof The proof relies on the adjngate identity adj(B) = det(B)B~^ applied to B = I— W. 
Recall that for a matrix B = ( 6 ij)ij=i,,,.,Ar, the adjngate of B is dehned as 

adj(B) = (det(B(^*)))^^^.^^_ (9) 

where is the matrix obtained by setting bji = 1, b^i = 0 for k j and bjk = 0 for A: 7 ^ i in B. 
Combining Eqnation Q with the identity adj(I—W) = det(I—W) (I—W)~^ gives for |||W||| < 1 

N i 

adj(I-W) = X Ew‘ = EE ( 10 ) 

k =0 k >0 e >0 k =0 


with fik = 0 for k > N. It remains to show that M = adj(I—W). We proceed entry-wise, hrst 
considering the case i = j. By constrnction of the adjngate matrix in ([^ for B = I —W, the 
conditions bu = 1 and bki = = 0 for k i correspond to setting Uki = tOik = 0 for all k. 

Plngging these valnes into det(I—W) = X)ceC 1° every cycle crossing Vi, yielding 


(adj(I-W)),, = //(c)c= ^//ii(/i)/i. (11) 

cec h&n 

i^V(c) 


Now consider the case i j. 
adj(I—W) satishes 


Going back to Eqnation (10), we see that the {i,j) entry of 


N 

(adj(I-W))^^. = Y'^kX ^ l^ic)fij{'w)cw. 

k=0 k>0 cSC uiSVV (c,io)sCxW 
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The right-hand side is a sum over hikes of the form h = cw with c a self-avoiding closed hike and 
w a walk from Vi to vj. By ([^, the left-hand side is obtained by plugging the values Uji = — 1, 
^ki = ^jk = 0 for k i,j and uu = cojj = 1 into det(I —W) = /^(c)c. Since only hikes from 

Vi to Vj remain (by identihcation with the right-hand side), we deduce in view of (|^, 

(adj(I-W)).^. = - ^ ^{hujji)h = ^ fVij{h)h. (12) 

heHij h&H 

Thus, adj(I-W) = M. □ 


The restriction of M to hikes of length i satishes, by identifying the terms of equal 
degrees in Lemma [2.2| 

. (13) 

In particular, commutes with W for all ^ > 0. If the digraph contains few self-avoiding hikes 
of length i, the matrix may have many zero entries. The construction gives in this case a 
non-trivial sparse matrix in the commutant of W. This kind of problems has some applications in 
practice, for instance in the study of random processes. In [2], the authors investigate conditions 
under which the transition kernel of a finite state Markov chain observed at random times can be 
estimated consistently. They show that sparsity conditions, arising from particular state transi¬ 
tions known to be impossible, suffice to recover the transition kernel when it commutes with a 
certain matrix for which an estimator is available. The identifiability conditions in this model 
rely on the existence of a sparse matrix in the commutant. Similar problems and applications 
are studied in [13 nans]. 


From the definition of Hij, one verifies easily that := satisfies 


£-1 


hesi 


(14) 


k=0 


for all i,j = This shows that the matrix enumerating, up to the coefficient ( —1)”'^^^ 

the self-avoiding hikes of length i for all pairs of vertices can be obtained as a polynomial of 
W. Because fJ-ij{h) is trivially zero when i(h) = N, one recovers Cayley-Hamilton’s theorem by 
setting in Equation (13). The general case also gives a direct proof of the identity 

£-1 


^£ = -]T.^ktr (W^-^), 


(15) 


k=0 


that (14) yields 


for which different proofs can be found in [T^ and m- To prove it using Lemma |2.2[ observe 

N 

= E E (-ir‘'’c = 


tr 


E (-1) 


_-\\ri{c)c ^ 


cecn5f, 


cecn5^ 
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noticing that each self-avoiding closed hike c appears exactly I = i{c) times when snmming over 
i. Hence, Eqnation (15) follows directly by compnting the trace on both sides of the eqnality 


mW = 

k=0 


3 Incidence algebra on hikes 

Becanse C is stable by mnltiplication, it forms a monoid with the empty cycle 1 as identity 
element. A natnral partial order on C arises from division: d ^ C divides c, denoted by d\c, if 
there exists c' € C snch that c = dc'. The closed hikes ordered by division form a locally finite 
partially ordered set, or poset. 

Definition 3.1 The reduced incidence algebra of closed hike is the algebra of real valued functions 
on C, endowed with the Dirichlet convolution, defined for f,g:C^M. by 

/*5(c) = , cGC. 

d\c 

In this definition, the snm is taken over all divisors d € C of c, inclnding the trivial cycle 1 
and c itself. The rednced incidence algebra is isomorphic to the algebra of formal series, endowed 
with mnltiplication. Indeed, for two fnnctions f,g:C^M 

X ='^f*9{c)c. 

^ cgC ' ceC ' ceC 

It follows that the Dirichlet convointion is associative, commntative and distribntive over addi¬ 
tion. The fnnction b defined on C by <5(1) = 1 and b(c) = 0 for all c 7^: 1 is the identity element 
for this operation as we have, for any fnnction f on C, f*b = 6*f = f. We refer to |18| for a 
more comprehensive stndy on this snbject. 

Extending this strnctnre of open hikes is slightly more complicated. The set of hikes TL is not 
stable by mnltiplication and for this reason, defining a division over TL leads to some difHcnlties. 
However, hikes are stable by mnltiplication with a closed hike so that {TL, .) forms an act over 
(C,.). Thns, the division by a closed hike can be extended to open hikes: for h € TL, d divides h 
if d £ C and h = dh' for some h' £ TL. Similarly, the Dirichlet convointion can be extended to 
/ : C —)> M and g : M by 

/ * 9{h) = ^ , h£TL. 

d\h 

Here again, the snm is taken only over closed divisors of h. Remark that / and g cannot be 
permnted in this expression, nnless h £ C. Interesting combinatorial properties arise from the 
poset strnctnre of hikes, by considering each hike individnally in Eqnation (|^. The dnality 
between walks and self-avoiding hikes becomes apparent in the next theorem. We extend the 
definition of fij to TL by setting fij{h) = 0 whenever h is not connected. 



Theorem 3.2 For all i,j = 1, N, (j,ij = * fij. 


Proof. From Lemma 2.2 we know that whenever |||W||| < 1, 

M = det(I-W) (I-W)-^ 


From Q and (|^ and using that det(I —W) = XlceC obtain for all i,j 

^ Hij{h)h = ^m(c)c X ^ fij{h)h = ^ /r* /^^(/i) h. 
h^'H cec h&'H h&'H 

The result follows by identification. 




□ 


This theorem reveals a somewhat unexpected relation between the function pLij which only 
takes non zero values for self-avoiding hikes and fij whose support contains only walks. Taking a 
closer look, the result is not surprising if h is self-avoiding. Indeed, a self-avoiding hike h from Vi 
to Vj has exactly one divisor d € C such that h/d is a walk from Vi to Vj. Thus, the convolution 
/i * fij is calculated over only one non-zero element and the equality is easily verified in this 
case. The result is actually more interesting if h is not self-avoiding as it yields in this case the 
non-trivial identity 

yhGn\s, =0. 

d\h 

Clearly, the function /i is a key feature to understand the combinatorial properties of this 
poset. The fact that /^(l) = 1 7 ^ 0 makes it invertible through the Dirichlet convolution, and its 
inverse /3 : C —)• M is the unique function characterized hy fi*P = fd*iJ, = 6. A reversed relation, 
expressing fij in function of iHj can then be derived easily, noticing that 

fij = {l3*n)* fij = fij) = 13* Hij. 

This relation turns out to be particularly important for our purposes, as we show that f3 satisfies 
interesting properties. In particular, we establish in the next proposition an expression of (3{c) 
that involves the number of appearances of each edge and vertex in c. Let Tij{c) denote the 
multiplicity of ojij in c and rj(c) = Aj(c) the number of edges starting from Vi in c (counted 

with multiplicity). 


Theorem 3.3 The function ft satisfies for all c ^ C, 

N 


m = n 


Ti(c)l 


Til(c)! X ... X TiN{c)\ 


Proof. We will prove that the function ft as defined in the theorem is the inverse of /i through 
the Dirichlet convolution. The case c = 1 being trivial, we take c ^ 1. Since //(d) = 0 if d is not 
self-avoiding, we have 

T* P{c) = '^p{d) 13(^2) = 

d\c d\c 

deCn5 
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For d ^ 1 a self-avoiding divisor of c, denote by the permutation over V{d) associated to d, 
i.e. such that d = n., ev{d)^iad{i)- Clearly, Tij{cjd) = Tij{c) if Vi ^ V{d), while for Vi £ V{d), 
we have Ti{c/d) = Ti{c) — 1 and 


It follows 



Tij{c)-1 iij = ad{i) 
Tij{c) otherwise. 


N 


c =n 


Ti{c/d)\ 


d' Til (c/d)! X ... X TiN{c/d)\ 


n 

Vi&Vid) 


Ti(c) 


X /3(c). 


We get, including the case d = 1 



E n 

d\c Vi£V{d) 

d£CnS\{l} 


Tiad{i){c) \ 

Ti{c) ) 


Recall that for d G Cn5^, /i(d) = (—1)^ sgn(iTrf). Let £ = |R(c)| (the number of different vertices 
in c), the previous equality becomes, regrouping the divisors with equal lengths. 


+ ^ sgn(ad) Yl 

d\c k=l d\c Vi£V(d) 

decns^ 


Ro-d(i)(c) \ 

n{c) ) 


Now consider the matrix B(c) with entries Tij{c)/Ti{c) for Vi,Vj G R(c). By identifying each 
self-avoiding divisor d of c with its corresponding permutation ad, we recognize in the above 
expression the characteristic polynomial of B(c) taken at A = 1, 


l + ^(-l)^ sgn(t7rf) Yl = det(I-B(c)). 

k=l d\c ied ' 

decn5'= 


Since B(c) is a stochastic matrix, det(I—B(c)) = 0 which ends the proof. 


□ 


The coefficient /3(c) corresponds to the number of arrangements of the edges in c, regrouped 
by their initiating vertex. Indeed, the multinomial coefficient 


Tijc)'. 

rji(c)! X ... X TiAr(c)! 


counts the ways of ordering the edges initiating from vi in c, accounting for their multiplicity 
Tij{c). Considering all configurations for each vertex in c recovers the coefficient /3(c). So, /3 
enumerates the different ways to travel along a closed hike. 
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This result points out some interesting properties of /?, most of which are not straightforward 
from its initial definition as the inverse of fj,. The first immediate consequence is that /? is positive. 
Secondly, /3(c) is equal to one if c is a self-avoiding closed hike. This condition is sufficient but 
not necessary, as we have for instance /3(c^) = 1 as soon as /3(c) = 1. A third consequence is that 
/3 is non-decreasing with respect to multiplication, which can be stated formally as: Vci,C 2 G C, 
/ 3 (ciC 2 ) > max{/3(ci),/ 3 (c 2 )}. Finally, j3 is multiplicative over decompositions on disjoint closed 
hikes. Indeed, if c can be written as the product of say p > 2 mutually vertex-disjoint compo¬ 
nents Cl, ...,Cp G C, then /3(c) = /3(ci).../3(cp). This property is reminiscent of the multiplicity of 
arithmetic functions over coprime integers (see for instance HI). In this framework, two closed 
hikes Cl, C 2 can be considered coprime if they share no common vertex. The multiplicity property 
of /3 is then inherited from the multiplicity of its inverse p. 


Remark. The function p is the Mobius function (i.e. the inverse of the constant function equal to 
1) of the trace monoid of circuits described in [^. Circuits correspond to partially commutative 
versions of closed hike so that, in our fully commutative framework, circuits composed of the 
same edges are seen as the same object. Thus, /3(c) counts the number of circuits composed of 
the same edges as c. 


A different expression for /3 can be derived from the inverse relation in Lemma 2.2 writing 
in terms of the A: = 0,1, This result is given as a corollary. 


Corollary 3.4 // 111 W| 11 < 1, then for £ G N, 


l 

(16) 

fc =0 

where the eoefficients 4>o, are defined by 

00 = 1 , ^ {-lyfik^-Tpkp, k = (17) 

ki-\-...-\-kp=k 

Before proving the result, let us clarify that the 0fc’s are defined by taking the sum over all com¬ 
positions of k, that is, all positive tuples (fci,..., kp) such that ki + ... + kp = k, for all p = 1,..., k 
(two tuples composed of the same integers ki,...,kp but in different orders are to be counted 
twice). 


Proof. This holds if W is the null matrix with the convention \N^ = I in this case. For W 7 ^: 0, 


Equation (10) yields 

m = Em‘'’ = E ipk X Ew 


N 


£>0 k=0 k>0 

We use the formal series expansion 

1 1 


Ew' 

£>0 




k>0 


Ek=o 1 + Ef=i A ^ 
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By regrouping the terms of equal degree, we get 


T.toi’k 




/c^l k\~\~.. ,~\~kp — k 


k>0 


Hence, 


£>0 fc>0 fe>0 


EE 4 M 

£>0 k=0 


and the result follows by identihcation. 


( 18 ) 


□ 


Like the V’fc’s, the coefficients 4>k are homogenous polynomials of degree k in the ujij’s. While 
^pk only involves the self-avoiding closed hikes, (/>*, depends on all the closed hikes of length k on 


the digraph. The formal series inversion in Equation (18) actually corresponds to the inversion 


of the Dirichlet convolution when identifying each closed hike. This means in particnlar that the 
coefficient (/)fc can be expressed as 

4>k= Y 

ceC'' 

One can verify this formnla directly from the formal series mnltiplication 

I = Y^f<^ ^ ^ Y = Y^* 


fc >0 


A:>0 


ceC 


cec 


cec 


recovering exactly the formal series version of the eqnality 6 
0 and ([T^, we dednce a new expression of /3(c) for c ^ 1: 


/3 * /X. By combining Equations 


«=) = E E (-i)V(ci)...Mcp) 

P>1 Ci...Cp=C 


E ( 

Si ...Sp —c 


2^n'(si)+...n'(sp) 


( 20 ) 


setting n'(.} = n(.) -|- 1, where the hnal sum is taken over all p-tuples (si, ...,Sp) of non-empty 
self-avoiding closed hikes snch that si...Sp = c. This eqnality provides an expression of /3(c) in¬ 
volving the different decompositions of c into self-avoiding closed hikes. While this expression is 
presnmably less practical than the previons one, it induces nevertheless interesting conseqnences 
from a combinatorial point of view, which are discussed in Section 


We now come to onr hnal resnlt, which expresses the mnltiplicity of an open walk in terms 
of its decompositions into self-avoiding components. This resnlt will be illustrated on some 
examples in Section]^ 


Theorem 3.5 Let h be a non-empty hike from Vi to Vj, 




E 




si...SkP=h 


setting n'{.) = n(.) -|- 1, where the sum is taken over all self-avoiding deeompositions of h, i.e., 
all {k -|- l)-tuples (si,..., Sfc,p) G (C n 5 \ {1})^ x Sij with k > 0 such that si...SkP = h. 
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Proof. Plug the expression of /3 in Equation (20) into fij{h) = Ylid\h P{d)fJ‘ij{h/d). For i ^ j, the 
fact that fj,ij{h/d) = 0 for hjdfi Sij simplihes into 

E X j^^n'(si)+---+n'(sfe)+n'(p) 

d=si...Sf^\h s\...SfiP=h 

h/d&Sij 

for p = h/d., thus recovering the result. For i = j, we use that = p{c)l{vi ^ F(c)}, 


fii{h) = '^(3{d)idi 
d\h 


E 

d\h 


E 


d\h 
ViGV{h/d) 


The hrst term of the right-hand side is —/3 =t= p{h) which is zero for all h ^ 1. The result follows 
by using the expression of j3 given in Equation (20), similarly as in the case i ^ j. □ 


4 Examples 

In this section, the functions fij,ij,ij,fi and /x are computed on some examples. For ease of 
comprehension, we start with explicit simple cases before considering more general structures in 
the hnal examples. 

Example 1. Let us begin with the graph represented in Figurej^which contains only two disjoint 
simple cycles. 



Figure 4: Disjoint simple cycles. 


The hike covering the whole graph is the closed hike hi = a;i2aJ23‘^34^^4i<^56^^67<^75 (recall that 
the order is not important). We obtain directly fij{hi) = 0 (because hi is not connected) and 
hijihi) = 0 (because hi crosses every vertex) for all i,j = 1,...,7. Moreover, the dehnitions 
of p and (3 give in this case p.{hi) = (3{hi) = 1. To check the equalities p * fa = nu and 
P* hii = fa, the calculations are straightforward, since the only closed divisors of hi are hi, ci = 
wi 2 W 23 W 34 a; 4 i, C 2 = and the void cycle 1. We get for instance, 

h * fii{hi) = pihi)fii{l) + ix{c2)fii{ci) = 0 = ^ii(hi), 
using that /ii(ci) = 1, /ii(c 2 ) = 0 and /x(c 2 ) = —1. From /xii(ci) = 0, we also verify 
13 * /Uii(hi) = /3{hi)pii{l) + /3{ci)pii{c2) = 0 = /ii(hi). 
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Example 2. We now consider the graph given in Figure composed of two simple cycles sharing 
one vertex and the covering hike /12 = 



Figure 5: Two simple cycles with one common vertex. 


Since /i 2 is closed, fij{h 2 ) is null for all i 7 ^ j. Moreover, there are two ways to travel across /12 
starting from V 2 , depending on which side is visited hrst, and one way for every other vertex. We 
deduce f 22 {h 2 ) = 2 and /ii(/i 2 ) = 1 for i = 1,3,4, 5. Since /12 is not self-avoiding /r(/i 2 ) = 0 and 
Theorem 3.3 gives (3{h2) = 2. To check the formulas, we now consider all the decompositions 
of /i 2 into a product of closed hikes. The two non-trivial divisors of h 2 are ci = W 12 W 23 W 31 and 
C 2 = ^ 24 ^ 45 ^ 51 . We verify for instance. 


h * f22{h2) — h(l)/22(h2) + li(c2)/22(ci) S" f 22 ( 02 ) — 0 — p^22(h2) 
13 * pLii(h2) = l3(h2)pLii(l) + l3(ci)pLii(c2) = 1 = /ll(/l 2 ) 


Example 3. This example deals with the walk = a;i 2 Ci; 23 ca 34 a; 4 ia; 62 <^ 25 ‘^ 54 W 46 composed of two 
simple cycles sharing two vertices, represented in Figure The non-trivial divisors of /13 are 
detailed in Figure 



Figure 6 : Two simple cycles with two common vertices. 


Clearly, fijih^) = 0 for i 7 ^ j. We compute the values fuihs) by enumerating all ways to travel 
across hs starting from Vi. For instance, the two ways from vi are 


1 ^ 2 ^ 3 ^ 4 ^ 6 ^ 2 ^ 5 ^ 4^1 

1 ^ 2 ^ 5 ^ 4 ^ 6 ^ 2 ^ 3 ^ 4^1 


We hnd fii{h^) = fsaihs) = hbihs) = fm(h^) = 2 and 722 (^ 3 ) = fii(h^) = 4. Here again, 
IJ-ijihs) is null for all i,j = 1,... ,6 as well as ^(h^), while ^(hs) = 4. We verify the convolution 
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equalities nuihs) = ^ * and fuihs) = (3 * naihs) for arbitrary vertices, e.g. 

fJ- * fiiihs) = /i(l)/ii(/i3) + A*(c2)/ii(ci) + At(c3)/ii(c4) = 0 = fiiiihs) 
f3 * ^i22{h‘i) = /3(/l3)/^22(l) = 4 = f22{h‘i) 



Figure 7: Simple divisors of / 13 : ci = a;i 2 a; 23 a; 34 a; 4 i (top-left), C 2 = a; 25 W 54 a; 46 W 62 (top-right), C 3 = 
wi 2 W 25 W 54 a; 4 i (bottom-left) and C 4 = a;23a;34a;46a;62 (bottom-right). 

Example 4- We consider the closed hike /14 = a;i 2 W 23 W 35 CU 56 u; 64 W 4 iW 25 W 54 a; 42 , illustrated in Figure 
composed of 2 cycles sharing 3 vertices. In this case, note that the orientation of the two cycles 
has an impact on the values of fa. 



Figure 8 : Two simple cycles with three common vertices. 
The walk /i 4 has 8 non-trivial divisors, detailed in Figure 
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Figure 9: Simple divisors of / 14 . From top-left to bottom-right: Ci = UJ 25 UJ 54 UJ 42 , C 2 = 

0Ji20J23<jJ35<^56<^64<jJ41, C 3 = U;i2W25W54W4i, C 4 = 0 J230J330J5e0Jg40J42, C 5 = a;25a;56a;64W42, Cg = 

W 12 W 23 W 35 W 54 W 41 , C 7 = ^23^^35^54^42 and Cg = a;i2W25W56W64W41- 


We know that fij{h 4 ) = 0 for i / j, iJ,ij{h 4 ) = 0 for all i,j = 1,... 6, ^(/i 4 ) = 0 and /3(/i4) = 8. 
Counting the connected paths starting from each vertex gives /ii(/i 4 ) = = fgg{h 4 ) = 4 

and 722 (^ 4 ) = / 44 (h 4 ) = = 8. We recover the correct values from the identities /4 m(/i4) = 

fJ' * fii{h 4 ) and fn{h 4 ) = (5 * Hii{h 4 ). Keeping only the non-zero values in the convolution gives, 
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for instance 


^22(^14) = /^(l)/ 22 (/l 4 ) + At(c 2 )/ 22 (ci) + f22(02) + ^(04) f22(03) + f^( 03 ) f22(04) 

+ /U(C 6 )/ 22 (C 5 ) + /U(c 5 )/ 22 (c 6 ) + /U(c8)/22 (c?) + /U(c 7 )/ 22 (c 8 ) 

=8-l-l-l-l-l-l-l-l=0 

/ll(/l4) = /3(/l4)/Ull(l) + /3(c2)/ill(ci) + /3(c3)/iii(c4) + /3(c6)/Uii (C5) + /3(c8)/Uii (C7) 

= 8- l- l- l- l = 4. 

Example 5. Consider a closed hike /15 composed of two simple cycles of opposite directions 
crossing n times. This example can be represented as n cycles placed one after the other. As 
we observed in the previous examples, the length of these cycles does not impact the values of 
the functions fa, pLu, fj, and /3 so that we can take cycles of length 2 without loss of generality, 
considering for instance the closed hike a;i 2 W 2 iW 23 a ;32 • • ■ uJni^in illustrated in Figure [T^ 



We find fij(h 3 ) = 0 for i / j, fn(h 3 ) = 2 n, ptij(h 3 ) = 0 for all i,j = l,...n, ^(h^) = 0 and 
/?(^5) = 2 ”'. Let ak = Ukk+i^k+ik, k = 1,... ,n — 1, an = iVni^in be the length 2 divisors of /15 
and do = uJi2i^23 • • • i^n-in^ni, di = 0^21(^32 ''' ‘^nn-ii^in the outside and inside cycles composing 
/15. The non-trivial divisors of /15 are do, di and every product Ofcj • • • Ukp obtained for a subset 
{ki,--- ,kp} of {l,...,n}. Using that /n is non-zero only for walks passing through vi and p, 
vanishes for non self-avoiding closed hikes, we obtain by keeping only the non-zero terms in the 
Dirichlet convolution 

n 

PL* fu(h5) = pi(l)fu(h) + pi(di)fii(do) + iJ-(do)fii(di) + '^Kak)fii 

k=2 

which recovers ultimately /x * /ii(/i 5 ) = 2n — 2 — 2(n — 1) = 0 = /xii(/i 5 ). The reverse relation 
/ii(h 5 ) = (5 * /iii(/i 5 ) is less trivial. To compute it, we have to enumerate for any p, the sets 
{ki, • • • , kp} C {1,..., n} such that /xii(afcj • • • 0 *,^) / 0 , i.e., such that is a self-avoiding 

closed hike that does not cross r;i. For each such closed hike, the complement h^Kak^ ■ • ■ akp) 
is composed of p disjoint connected components, one of which contains aia^. This component 
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can be divided into two connected components by separating oi and Thus, each closed hike 
Ofci • • • such that • • • akp) 7 ^ 0 can be associated with a composition of n —p containing 

p + 1 elements, which there are of them. For each • • - akp, the coefficient /3 of the 

complement equals 2 ”“^^ and one recovers the formula 

[ "-i j 

P*m{h5)= iz 2"“''’ = 2n = /ii(h5). 

p=0 \ P / 


Example 6. Consider a self-avoiding closed hike composed of n > 2 simple connected compo¬ 
nents ai, ...,an (we may assume without loss of generality that each connected component is a 
loop as illustrated in Figure 11). 



From the first expression of (3 given in Theorem 3.3 it is clear that = 1. On the other hand, 

the number of ways to decompose hg into a product oi k < n non-empty self-avoiding closed 
hikes writes as the sum of the multinomial coefficients over all positive compositions ni,..., of 
n. Since for any self-avoiding decomposition si,...,Sk, the product p(si)...p(sfc) always equals 
fJ-ihe) = (— 1 )”, combining the two expressions of (3[ho) yields 


/ 3 (/! 8 ) = 1 = {-!)” (- 1 )‘ 


n! 


k=l 


n-i+...+nk=n ^ 


Alternatively, this equality can be obtained by identifying the coefficient of in the power series 
expansions of e~^ = 1/e^. Since /ig is not connected, we know that fuihg) = 0 for all i. To 


compute the expression of /ii(hg) from Theorem 3.5 we consider the self-avoiding decompositions 
of the form hg = si...SkP with si,..., Sk € C CiS \ {1} and p G Su. To verify that this expression 
gives /n(/ig) = 0 in this case simply observe that any self-avoiding decomposition si,Sk,p 
such that p 7 ^ ojii cancels out with the decomposition si,..., Sk,pjoJii,u]ii in view of 


j'_]_^n'(si)+...+n'(sj,)+n'(ui) _ {si)+...+n'(s^+n'{p/LJii)+n'(uju) 
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Thus, summing over all self-avoiding decompositions recovers /^(hg) = 0. 


Example 7. We now consider the closed hike hj constructed from the previous example with an 


where cq = a;i 2 W 23 <^ 34 W 45 W 56 a; 67 a; 78 a; 8 i)- In this example, the cycle cq is isolated in every self¬ 
avoiding decomposition since it shares a common node with all the other divisors of / 17 . Thus, 
the different ways to express h-j as a product of self-avoiding closed hikes can be obtained from 
the previous example, inserting the cycle cq wherever possible. Precisely, for a decomposition 
/ig = si...Sfc of hg into k < n non-empty self-avoiding closed hikes, there are exactly k + 1 
possibilities to insert cq. Moreover, remark that fi{si)...fi{sk)fi{co) = /i(/ig)/r(co) = (—is 
constant over all self-avoiding decompositions. Combining the two expressions of j3{hj) thus 
recovers the formula 

n 

(_i)‘+i(fc+i) ^ (_ 

k=l ni+...-|-nfe=n 



1 ) 


n+1 


n\ 


ni!...nfc! 


= 2^ 


extra cycle cg passing through each vertex: hj = Hq x cg (e.g. in the right graph in Figure 11 


In this example, there are two ways of visiting the whole walk hj from one vertex Vi to itself, 
depending on whether the loop at Vi is traveled at the start or at the end. Thus, /n(/i 7 ) = 2 for 
all Vi- In a self-avoiding decomposition with si,..., G Cn5\{l} and p G Su we can distinguish 
the cases p = cq, p = uju and p 7 ^ tCjj, cg. Clearly, the sum over all self-avoiding decompositions 
/17 = si...,Sfcp such that p = cg yields /3(hg) since (—1)” = 1. Moreover, the sum over all 

self-avoiding decompositions with p = oJa recovers jujii) = 2”'“^. Finally, for a self-avoiding 

decomposition hg = s\...S]^p of /ig with p ^ eoa, there are k possibilities to insert cg, yielding 

= /3{he) + /3(—)+ £(-1)"=+'^: = 1 + 2-1 - 2-i + 1 = 2. 


Acknowledgments 

The authors are grateful to Pierre-Louis Giscard for his explanations on the poset structure 
of hikes, and to an anonymous referee for its helpful comments which helped improve this paper. 


References 

[ 1 ] Tom M Apostol. Modular functions and Dirichlet series in number theory. AMC, 10:12, 
1990. 

[2] Flavia Barsotti, Yohann De Castro, Thibault Espinasse, and Paul Rochet. Estimating the 
transition matrix of a Markov chain observed at random times. Statistics & Probability 
Letters, 94:98-105, 2014. 

[3] Philippe Blanchard and Dimitri Volchenkov. Random Walks and Diffusions on Graphs and 
Databases: An Introduction, volume 10. Springer Science & Business Media, 2011. 


19 





[4] Raffaella Burioni and Davide Cassi. Random walks on graphs: ideas, techniques and results. 
Journal of Physics A: Mathematical and General, 38(8):R45, 2005. 

[5] Peter J Carrington, John Scott, and Stanley Wasserman. Models and methods in social 
network analysis, volume 28. Cambridge university press, 2005. 

[6] Pierre Cartier and Dominique Foata. Problemes combinatoires de commutation et rear¬ 
rangements. Lecture notes in mathematics, 85, 1969. 

[7] Fan RK Chung. Spectral graph theory, volume 92. American Mathematical Soc., 1997. 

[8] Dragos M Cvetkovic. Eigenspaces of graphs. Number 66. Cambridge University Press, 1997. 

[9] Xingguang Han, Yangzhou Chen, Jianjun Shi, and Zhonghe He. An extended cell transmis¬ 
sion model based on digraph for urban traffic road network. In Intelligent Transportation 
Systems (ITSC), 2012 15th International IEEE Conference on, pages 558-563. IEEE, 2012. 

[10] Frank Harary. The determinant of the adjacency matrix of a graph. Siam Review, 4(3):202- 
210, 1962. 

[11] Zhang Hong-Hao, Yan Wen-Bin, and Li Xue-Song. Trace formulae of characteristic poly¬ 
nomial and cayley-hamilton’s theorem, and applications to chiral perturbation theory and 
general relativity. Communications in Theoretical Physics, 49(4) :801, 2008. 

[12] Yasunari Inamura. Estimating continuous time transition matrices from discretely observed 
data. Bank of Japan Working Paper Series, (06), 2006. 

[13] Gregory F Lawler. Loop-erased self-avoiding random walk and the laplacian random walk. 
Journal of Physies A: Mathematieal and General, 20(13):4565, 1987. 

[14] Elizabeth Chase MacRae. Estimation of time-varying Markov processes with aggregate 
data. Eeonometrica, 45(1):183-198, 1977. 

[15] John S Maybee, DD Olesky, Driessche P van den, and G Wiener. Matrices, digraphs, and 
determinants. SIAM Journal on Matrix Analysis and Applications, 10(4):500-519, 1989. 

[16] A. O. Pittenger. Time changes of Markov chains. Stochastic Process. AppL, 13(2):189-199, 
1982. 

[17] J Ponstein. Self-avoiding paths and the adjacency matrix of a graph. SIAM Journal on 
Applied Mathematics, 14(3):600-609, 1966. 

[18] Gian-Carlo Rota. On the foundations of combinatorial theory. In Classic Papers in Combi¬ 
natorics, pages 332-360. Springer, 1987. 

[19] Loth Asker Zadeh and Charles A Deoser. Linear system theory. Robert E. Krieger Publishing 
Company, 1976. 


20 



